arXiv:1506.07887vl [astro-ph.IM] 25Jun2015 


Astronomy & Astrophysics manuscript no. FP_arxiv ©ESO 2015 

June 29, 2015 


Calibrating echelle spectrographs with Fabry-Perot etalons 

F. F. Bauer, M. Zechmeister and A. Reiners 

Institut fur Astrophysik (IAG), Georg-August-Universitat Gottingen, Friedrich-Hund-Platz 1, D-37077 Gottingen 
e-mail: fbauer@astro. physik. uni-goettingen. de 

submitted 


ABSTRACT 

Context. Over the past decades hollow-cathode lamps have been calibration standards for spectroscopic measurements. Advancing 
to cm/s radial velocity precisions with the next generation of instruments requires more suitable calibration sources with more lines 
and less dynamic range problems. Fabry-Perot interferometers provide a regular and dense grid of lines and homogeneous amplitudes 
making them good candidates for next generation calibrators. 

Aims. We investigate the usefulness of Fabry-Perot etalons in wavelength calibration, present an algorithm to incorporate the etalon 
spectrum in the wavelength solution and examine potential problems. 

Methods. The quasi periodic pattern of Fabry-Perot lines is used along with a hollow-cathode lamp to anchor the numerous spectral 
features on an absolute scale. We test our method with the HARPS spectrograph and compare our wavelength solution to the one 
derived from a laser frequency comb. 

Results. The combined hollow-cathode lamp/etalon calibration overcomes large distortion (50 m/s) in the wavelength solution of the 
HARPS data reduction software. Direct comparison to the laser frequency comb bears differences of only maximum 10 m/s. 
Conclusions. Combining hollow-cathode lamps with Fabry-Perot Interferometers can lead to substantial improvements in the wave¬ 
length calibration of echelle spectrographs. Etalons can provide economical alternatives to the laser frequency comb, especially for 
smaller projects. 

Key words. Instrumentation: interferometers - Instrumentation: spectrographs - Methods: data analysis - Techniques: radial veloc¬ 
ities 


1. Introduction 


Accurate wavelength calibration is a cornerstone for any mea¬ 
surement with high resolution spectrographs. It translates the 
pixel position of the detector to absolute wavelengths and 
thereby defines the physical scale. As many astronomical 
research fields benefit from high precision radial velocity 
(RV) measurements, the next generation of spectrographs aim s 
for cm/s precisions (iD’Odorico & CODEX/ESPRESSO Teaml 
2007h . Today, most echelle spectrographs use hollow-cathode 
lamps (HCLs), e.g., Thorium-Argon or Thorium-Neon, to per¬ 
form wavelength calibration but these lamps are insufficient for 
the needs of future high precision instruments. 

Hollow-cathode lamps establish collisions between the car¬ 
rier gas (Argon or Neon) and the Thorium atoms which thereby 
get exited (iKerber et al.1120071) . When the Thorium relaxes into a 
lower energy state, a photon carries away the energy difference 
and contributes to an emission line. The laboratory wavelengths 
of Thorium lines can be measured or ca l culate d as done by 
Palmer & Engl emanl ( 19831) . Lovis & Pend (120071) . IKerber et~al] 
d2008l) and iRedman et al. (12014 ). The drawback of HCLs is 
the limited number of reference lines, the large dynamic range 
spanned by the lines, and the irregular line distribution. Blend¬ 
ing of lines is an additional problem: in total only 4000 - 8000 
spectral lines remain usable in the visual dLovis & Pepeil2007 ). 
In addition to the large intensity differences between Thorium 
lines, the carrier gas often emits strong lines itself. Exposure 
times have to be well chosen in order to have enough cathode 
(e.g. Thorium) lines for calibration while keeping the number 
of saturated carrier gas (e.g. Argon) lines low. Spectral regions 


without lines must be interpolated by the wavelength solution 
resulting in uncertain wavelengths for these pixels. 

To overcome the problems of HCLs new calibrators are 
currently under develop ment. Among them are the laser fre¬ 
quenc y comb (LFC, e.g. lMurphv et al.ll2007l [wilken et al.|2010. 
20121) and the Fabry-Perot i nterferometer ( FPL Wildi et al.l201C ; 
Schaf er & Reinersf 120121 : iReiners et alJ l2014t ISchwab et a 1 


20141) . LFCs are rather expensive and often not affordable for 


smaller projects. Here we focus on less cost intensive FPIs and 
their use as a calibration source. 


A tunable FPI device was already used as dispersing element 
to develop a high resolution spectrograph for planet search lo¬ 
cated at Kitt Peak dMcMillan et al.fl988l.ll994l) . The wavelength 
solution of the spectrograph relied totally on geometrical cali¬ 
bration and accurate tuning of the FPI demonstrating precisions 
of 6 m/s. [Cochran et alJ (119821) used the FPI to superimpose 
absorption lines of known wavelength in the stellar spectrum. 
Both used the FPI in transmission resulting in very low instru¬ 
ment efficiency leading to the development of new techniques 
such as iodine cells in tran smissi on (IButler et ali 1 1996 ) or ul¬ 
tra stable spectrographs (iMavor et all2003l) '. Currently, FPIs are 
again considered in astronomical spectro graphs for nightly drift 
checks, but not as wavelength calibrators (IWildi et al.ll2010h . 

The advantage of the FPI against HCLs is its ability to pro¬ 
duce a dense grid of lines with almost uniform intensity over 
the entire spectral range. They can be manufactured to match the 
specifications of most spectrographs as the free spectral range 
is a function of the cavity width. This tackles the problem of 
empty spectral regions and saturated lines thereby increasing the 
robustness of the wavelength solution. 
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wavelength A. [nm] 
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Fig. 1. Spectrum of the CARMENES NIR FPI taken with the Fourier 
Transform Spectrometer in Gottingen. 

The drawback of FPIs is their weakly constrained absolute 
wavelength. Line positions are determined by the cavity width of 
the etalon and the interference order. Both are usually not known 
to the accuracy needed to provide absolute calibration. Further¬ 
more, the group-velocity dispersion limits the predictability of 
individual peak positions, and the cavity width is very sensitive 
to temperature and mechanical stress. Hence FPIs_are tempera- 
ture and pressure stabiliz ed to minimize drifts (IWildi et al.fe o 10; 
iSchafer & Reinersll2012l) . Such a stabilized FPI is located at the 
HARPS spectrograph and nightly drifts of 10 cm/s were reported 
(tWildi et al.ll201Ch . 

First, we propose a method to characterize the FPI comb pat¬ 
tern with HCLs, using a high resolution spectrograph in Sect 2. 
We explain our wavelength calibration concept and model in 
Sect. 3. In Sect. 4 we combine information from an HCL and 
a FPI to improve the accuracy and stability of the HARPS wave¬ 
length solution as compared to the use of an HCL alone. Fur¬ 
thermore, we investigate the impact of a distorted wavelength 
solution on RV measurements in Sect. 5 

2. Calibrating the Fabry-Perot interferometer 

2.1. The ideal FPI 


the effective FPI cavity width d. We use the interference condi¬ 
tion of the FPI, Eq. ([[}, to perform the cavity width measure¬ 
ment. This can be done with any high resolution s pectrograph 
and its internal calibration sources such as HCLs (IWildi et al.l 
l20lol) . First, we derive the wavelength solution of the spectro¬ 
graph using only HCL lines. Then we derive the line positions 
of the FPI on the detector and use the HCL wavelength solu¬ 
tion to assign the corresponding wavelength to the interference 
peaks. 

Now we need to identify the interference order m for all FPI 
peaks. From the observed spectra we can only obtain a relative 
numbering k — m - m\ by simply counting the lines but the 
absolute interference order of the first (reddest) peak m\ is not 
known. To find m\ we use the fact that the ideal FPI cavity width 
is a constant number. We start with a guess for m\ and use the 
resulting interference orders m = k + m\ and the wavelengths de¬ 
rived for each line in Eq. ([TJ to obtain the effective cavity width d 
for each FPI line separately. If we plot the effective cavity width 
d as a function of relative peak number k using the correct m\ we 
expect the data points to be distributed around a constant value. 
If our guess for ni\ is wrong the data points will, however, pro¬ 
duce a positive slope for m\ too small or a negative slope for m \ 
too high. Once the right value for m\ is found the large number 
of individual effective cavity width measurements (typically in 
the order of 10 4 for high resolution echelle spectrographs) can 
be averaged to obtain a good estimate of d. 

We now only need this one parameter (the effective cavity 
width d of the FPI) to globally predict the wavelength of all inter¬ 
ference peaks observed. The combined accuracy of all HCL lines 
anchors the FPI on an absolute wavelength scale. Paired with the 
high precision of the FPI local imperfections of the wavelength 
solution can now be resolved and corrected. The numerous FPI 
features bridge regions lacking any spectral features of HCLs 
and we can now check if the wavelength solution model is suffi¬ 
cient or if more detailed models are required. 

2.2. The real FPI 


The ideal FPI consists of two partly reflective plane parallel sur¬ 
faces dBorn & Wolf! 11999 ). When light enters the device it gets 
reflected between the two surfaces and interference takes place 
between light rays depending on the phase difference. With a 
white light input source, a spectrograph will observe a pattern of 
transmission peaks spaced equidistantly in frequency v (Fig.Q}. 
The transmission maxima fulfill the interference condition given 
by the effective FPI cavity width d: 

mA m = 2d (1) 

where m is an integer number called interference order and A m 
are the transmission maxima wavelengths. The effective cavity 
width is d = Dn cos 6 with the distance between the two mirrors 
D, the refraction index n and the incident ray angle 6. 

Using the FPI as a calibration source requires knowledge of 
the exact wavelengths for all interference maxima observed by 
the spectrograph. As the wavelength A n , of each FPI interference 
order depends only on the effective cavity width, d, this param¬ 
eter must be known with the same accuracy as the wavelength 
requirement (6d/d = 3 • 10 9 for 1 m/s). For the FPI this means 
that the spacing d must be known with an accuracy of about 
0.15 A dReiners et al.ll2014)i . Typically the mirror distance D is 
only known to about 1 pm which is five orders of magnitudes 
above our requirements. In addition we do not exactly know the 
incident ray angle 8. Hence, we need to find a method to calibrate 


In practice, applying the concept of calibrating the effective cav¬ 
ity width using standard calibrators proves to be more difficult 
for real FPIs. The cavity width is not constant throughout the 
wavelength range covered by high resolution echelle spectro¬ 
graphs. FPIs are usually soft coated allowing photons of differ¬ 
ent energy to penetrate to different depth of the dielectric sur¬ 
face. Hence, photons of different wavelength see different cav¬ 
ity widths. As an example the penetration depth variation over 
the wavelen gth ra nge of HARPS easily reaches hundred nm 
(IWildi et al.l 2010 ). If we do not account for this the computed 
FPI wavelengths are incorrect by several km/s. For a full charac¬ 
terization of the FPI spectrum, the task is therefore not only to 
find the global mean effective cavity width d of the FPI but also 
to determine the wavelength dependent function d(A). 

We assume that the penetration depth is a smooth function 
of wavelength and therefore also a smooth function of relative 
interference order k. Because the measured wavelength is uncer¬ 
tain but the relative peak numbering is noise-free, we model d 
as a function of k rather than d as a function of A. W e choose 
uniform cubic B-splines (Ide Booiil200ll : lDlerckxlll995h to fit the 
function d(k). 

The assumption that the penetration depth variation is 
smooth is the cornerstone of our method. Hence we decided to 
test this with the CARMENES ( Quirrenbach et al .l201 111 near in¬ 
frared FPI (finesse ~ 8. iSchafer & Reinersll2012 ). We obtained 
78 spectra with the Bruker IFS 125/HR Fourier transform spec- 
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Fig. 2. Upper panel: Effective cavity width measurement d{k) for the 
CARMENES NIR FPI. Lower panel: Residuals of spline fit. Data (black 
dots), clipped lines due to atmosphere absorption features (red crosses), 
spline fit to the data (green line). 


trometer (FTS) with a resolution of R = 500 000 in our optics 
laboratory. The median signal-to-noise ratio of the coadded spec¬ 
trum is 350 (Fig. |T1). As starting guess for all FPI peak wave¬ 
lengths we used the internal wavelength calibration of the FTS. 
It is provided by a He-laser measuring the path length within the 
spectrometer. With the high resolution of the FTS the FPI peaks 
are well resolved in the spectra and we measured individual line 
positions by fitting a Lorentz function. The FPI itself was not op¬ 
erated under vacuum conditions, thus water lines contaminated 
the spectrum. We excluded contaminated peaks from our cavity 
width measurement. We determined the mean effective cavity 
width d = 12.338 350 ± 0.000040 mm which is consistent with 
what we expect from the mirror distance given by the manufac¬ 
turer D = 12.334 mm, the refractive index of air n ~ 1.003, and 
an incident ray angle of 6 — 0 deg. 

The penetration depth variation of the FPI as a function of 
the relative peak numbering and the measured FTS wavelength 
is shown in Fig. [2] Fitting a B-spline with 50 nodes results in an 
RMS scatter around the model of 0.326 nm translating to 8 m/s. 
We could not reveal any discontinuities in d(k) within our mea¬ 
surement precision. Hence our assumption that the penetration 
depth variation is smooth is valid. The RMS scatter of 8 m/s is 
about what we expect from photon noise. 

Once we have determined the model parameters for the ef¬ 
fective cavity width d we need to find the absolute interference 
number of first observed FPI peak m\ to compute the wavelength 
for all interference orders. This is no longer trivial if the effective 
cavity width d is not a constant. 

iReiners et al.l (120141) suggested to use the mean free spectral 
range (Av) to calculate the mean effective cavity width and then 
use Eq . {!]) to calculate the absolute order numbering. [Wildi et al.l 
(1201 Of) determined the absolute order numbering assuming the 
effective cavity width d to be exactly the absolute cavity width 
D specified by the manufacturer and rounding m to the closest in¬ 
teger number they obtain from equation Eq. CQ> and the HARPS 
wavelength solution. Their approach resulted in the penetration 
depth curve with the minimum variation. 

Other choices of in | introduce steeper slopes in the effective 
cavity width d. We show how the choice of m \ influences the 
overall slope of the function d(k) for the CARMNENES NIR FPI 
in Fig. [3] Because we do not know enough about the material 
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Fig. 3. Upper panel: Effective cavity width d(k) obtained for different 
absolute interference order numbers ni\. The reference curve (mi = 
13604, black) is shown in Fig. [2] Lower panel: Difference in modeled 
FPI peak wavelength A m using m i = 13704 with the reference curve 
(dashed line). Residuals of the spline fit to the cavity width d{k) ob¬ 
tained with = 13704 (purple dots). 


properties of the coating, the penetration depth variation with 
wavelength and the alignment of the FPI, we do not know which 
m\ and which d curve are the true ones. We chose m\ = 13604 
as reference value because it results in the d{k) curve with the 
minimal penetration depth variation. We note that, this choice is 
not necessarily the true value of m \ 

We realized, however, it is not necessary to know m\ ex¬ 
actly to calibrate the FPI wavelengths. Due to the degeneracy 
in Eq. £[} between the absolute peak number m and the effective 
cavity width d(k), identical peak wavelengths occur for differ¬ 
ent combinations of m and d(k). Hence, our method to derive 
absolute wavelengths for individual FPI peaks should be robust 
against choosing the wrong value for mi. 

To test this we computed d(k) using two different values 
for mi differing by 100 interference orders, mi = 13604 and 
mi = 13704. For both cases we obtain the FPI wavelengths us¬ 
ing Eq. (jj}. We calculate the difference of the calibrated FPI 
wavelengths line by line and show the result as dashed line in 
the lower panel of Fig. [3 The maximum difference between the 
calibrated peak wavelengths is 10 cm/s in the region between 
k = 2000 and k = 10000. At the edges the difference increases 
up to 5 m/s due to the decreasing number of data points to con¬ 
strain the spline fit. Despite the large difference in d(k) (several 
hundred nm, see upper panel in Fig. [3 the difference between 
the modeled wavelengths is small (lower panel in Fig. [3- 


3. Wavelength solution model 

In this section, we explain our strategy to determine the wave¬ 
length solution of a high-resolution echelle spectrograph. First, 
we introduce our concept to carry out the mapping between 
wavelength and pixels in the case of a monolithic detector with 
uniform pixels. We motivate our model, parametrization and re¬ 
gression type. Second, we advance to real detectors and show 
how we account for uneven pixel sizes in the wavelength solu¬ 
tion. 
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3.1. Mapping wavelengths to detector coordinates 

The general idea of wavelength calibration is to feed the spec¬ 
trograph with a reference spectrum of known wavelengths A/, 
for example the known atomic lines of a hollow-cathode lamp, 
and observe the response at the detector. The positions of the 
line centers x/ in the extracted spectrum are measured by fitting 
a model of the instrumental profile (e.g. Gaussian functions) to 
each individual spectral line /. For the set of wavelength-pixels 
pairs (/!/, x/), a relation is established via a fit. Finally, this re¬ 
lation is used to predict the wavelength at the center of every 
detector pixel. 

The wavelengths of standard lines are given by accurate and 
precise line lists with uncertainties in the range of 10 6 - 10~ 7 
(~ 10 m/s) ( iRedman et al .i 2014 ). On the other hand, the posi¬ 
tions of weak lines on the detector exhibit typical measurement 
errors of 0.1 pixels or about lOOm/s. Hence, in wavelength cali¬ 
bration the standard line wavelength. A/, represents the indepen¬ 
dent variable (or the cause) and the position on the detector, x/, 
is the dependent variable (or the spectrograph response). 

A common and simple procedure to do wavelength calibra¬ 
tion is to fit a continuous model for A(x). From the function A(x) 
one can evaluate directly the wavelength for each pixel. In this 
approach the role between dependent and independent variable 
is reversed (llsobe et al.l 1990t) and we call this reverse regression. 
Reverse regression causes problems when weighting the data. 

For this reason one could use direct regression and fit the 
detector position of individual lines as a function of wavelength, 
x(A). To predict the wavelength for each pixel in this approach 
the model x(A) must be inverted A(x) = x(A)~ . 

The wavelength solution is usually modeled with polynomi¬ 
als which, in case of reverse regression is oA = poly(x) (e.g. 
iBaranne et al.l dl996l) l. can be interpreted as a Taylor expansion 
of the grating equation 

oA = ncr [cos y\ sin a\ + cos 72 sin a 2 ] (2) 

in terms of the diffraction angle 02 which is proportional to the 
position in the focal plane, x oc 0 : 2 - In Eq. 0, o is the diffraction 
order, n is the refractive index, <r is the grating constant, a 1 and 
a 2 are the angles of the incident and diffracted ray with respect to 
the grating normal, and yi and yi are the off-plane angles before 
and after the grating, resp. dSchroederll2000h . 

In direct regression with x(A) = poly(o/l) we basically ap¬ 
proximate a re-arranged version of the grating equation 


/ 1 1 

r °A ■ 1 \ 

ci2 = arc sin - 

-cosyi sinofi 

v cos 72 

ncr 


We note that polynomials of higher degree (typically +1) are 
needed with direct regression for a wavelength solution with 
similar quality compared to reverse regression. This is because 
the sine function in Eq. 0 is more amenable to a Taylor expan¬ 
sion than the arcsine function in Eq. (0. 

A way of keeping both advantages (lower number of param¬ 
eters and direct regression) would be to use x = poly(oAy 1 as 
forward model which we might call inverse regression. This re¬ 
quires non-linear least square fitting and thus higher computa¬ 
tional effort. 

In what follows we chose direct regression. This allows us to 
correctly weight all data points and to implement unequal pixel 
sizes in our model (next section). Moreover, instead of simple 
polynomials for each order, we use 2D polynomials to couple 
the individual orders. This decreases the number of parameter 
and increases the robustness of the fit. We also use oA instead of 


A as variable, i.e. x(A, o ) = poly(o/l, o). The re-parametrization is 
motivated by Eq. 0. This approach is a direct regression version 
of the algorithm impleme nted in the IDL REDUCE package of 
(tPiskunov & Valentifl2002l) . 


3.2. Unequal pixel sizes 


Standard wavelength calibration procedures usually assume that 
detector pixels are equally spaced and of the same size. Irregu¬ 
larities in pixel size of CCD detectors are on the order of 10 2 
(IWilken et al.l 12010 1. which is becoming a significant limitation 
for high-precision spectroscopic measurements. For example, 
IWilken et al.l (i2010i) showed that inhomogeneous pixel sizes lead 
to discontinuities in the wavelength solution of HARPS on the 
order of several 10 m/s. 

In order to account for irregular detector pixels, we include 
the pixel size in our model and we distinguish now more strictly 
between detector coordinates, X (pixel), and focal plane coordi¬ 
nates, x. We denote the size of pixel i by a plx ,- and we assume 
that there is no gap space between individual pixels. Thus, the 
transformation from pixel coordinates to focal plane coordinate 
is 


IX] 


X(X, flpjx) — Opix.i 6!pix,[X] ' I X | X \ 


(4) 


where [A] is the value of X rounded to the nearest integer (the 
pixel number). In Eq. 0 we integrate over the width of all 
preceding pixels (including the current one) and linearly inter¬ 
polate (backwards) for sub-pixel position. The zero point of 
the coordinate system (x = 0 ) coincides with the left border 
of the first pixel (X = 0.5). If all pixels have the same size 
(a p ix,/ = 1) the transformation is simply x = X - 0.5. The 
transformation can also be extended to include virtual pixels for 
the application to mosaic detectors with gap sizes correspond¬ 
ing to several pixels (fl n ; x ,■ ~ 100 pixels) between the detectors 
as, e.g., in CA RMENES (Ou irrenbach et al.1l201 1 1) or CRIRES+ 
(iFollert et al,ll20l4 . 

Ideally, the sizes of all pixels are known prior to wave¬ 
length calibration, but in practice the individual pixel sizes are 
not known precisely enough. For future detectors, it might be 
interesting to determine pixel sizes with a microscope prior to 
integration in the instrument. 

Determining the sizes of all pixels is beyond the scope of 
ordinary wavelength calibration because the number of free pa¬ 
rameters exceeds the number of comb lines which are separated 
by a few pixels. But for a few, significantly deviating pixels it is 
indeed possible to simultaneously determine pixel sizes and the 
wavelength solution, if a spectral calibration source with dense 
line combs (FPI or the LFC) is available. 

For example the HARPS detector exhibits pixel irregularities 
every 512 pixels (IWilken et al.ll 2010 t> . and we therefore choose 
the following pixel size model 


Wi, if modulo(f, b ■ 512 - 3) = 0 
1 , otherwise 


(5) 


where every 512th pixel is allowed to depart from unity, while 
all others are assumed to be homogeneous and are fixed to unity. 
The first pixel size irregularity, b = 1, is located in pixel 509 
motivating an offset of 3 pixels in Eq. 0. One can further limit 
the number of free parameters when using column sizes instead 
of individual pixels as we do in Sect. 4 for HARPS. 

To find the pixel sizes w,, we perform a non¬ 
linear least-square fit by minimizing x 2 = E/“ 2 K ~ 
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X(x(A/, oi, a potyj, flpi x )] 2 . To evaluate the x~ with the observed 
data Xi we must invert Eq. ©. The back-transformation of the 
focal plane coordinates x into the detector coordinates X is 


X(x, «pix) 


* - X j 

a pix,j 



( 6 ) 


where j is the pixel number into which x falls which can be found 
via the condition Xj < x < Xj+\ where Xj = x(j - 0.5,a P j X ) = 

2/ = l a pix,i • 

The best fitting model and coordinate transformation is then 
inverted numerically to obtain the wavelength for the center of 
each pixel. We found that a simple bisection method is appropri¬ 
ate 0. 

To accelerate the non-linear least-square fit we separate it 
into a linear and a non-linear least-square fit. For a trial set of the 
non-linear parameters a p i X , the position of all calibration lines 
can be transformed to focal plane coordinates x/(Xi ) with Eq. ((4}. 
Then a smooth 2D polynomial is fitted with a linear least-square 
fit. The smooth model is transformed back to detector coordi¬ 
nates with Eq. © and the x 2 with the observed data X\ is com¬ 
puted. Testing different trial sets (e.g. with a downhill simplex; 

I Press! (120071) f the best fitting pixel sizes a p i X are found. 


4. Calibrating HARPS using its FPI 

4.1. Detector characterization and wavelength solution 

To show the potential of Fabry-Perot interferometers in wave¬ 
length calibration procedures for echelle spectrographs, we cal¬ 
ibrate the HARPS spectrograph located at the ESO 3.6 m tele¬ 
scope in La Silla Chile dMavor et al.l12003!) with its etalon and 
our wavelength solution method. HARPS has a resolution of 
115 000 (one resolution element is sampled with about 3.2 pix¬ 
els), covers the wavelength range from 380 to 690 nm and offers 
ThAr, FPI a nd LFC spe ctra. The design of the HARPS FPI is 
described in IWildi et al.l (i2010i) . Briefly, the HARPS FPI has a 
cavity width of about 7.3 mm and a measured finesse of about 
4.3. It is fiber coupled and temperature and pressure stabilized 
inside a vacuum tank. 

We analyzed calibration data taken on March 16, 2011, of¬ 
fering both ThAr and FPI exposures. We reduced the calibration 
spectra using the IDL REDUCE package. The pixel position of 
the Thorium lines were measured by fitting a Gaussian profile to 
the extracted ID spectra. The FPI peaks are separated by 8 - 21 
pixels (from blue to red orders). The median FWHM of the FPI 
lines is about 5 pixels or 1.6 resolution elements which means 
the FPI lines are moderately resolved in the spectra. Hence the 
line shape of the FPI deviates from the Gaussian like instrumen¬ 
tal profile of HARPS and we fit a Lorentz function to each peak 
to obtain the pixel position. 

The first step in our calibration procedure is to derive the 
wavel ength solution from the Tho rium lamp only using the line 
list of iPalmer & Englemanl (119831) . Because of the low number 
of Thorium lines available in single calibration images we use a 
2D polynomial in the direct regression (7th degree in dispersion 
direction and 6th degree in cross dispersion) and uniform pixel 
sizes. The residuals of the best fit to the Thorium lines in block 
2 of the blue detector (echelle orders 116 to 135) are plotted in 
the upper panel of Fig. |4] 

2 Methods involving first derivatives are not recommended by us be¬ 
cause of convergence problems with discontinuous functions as it is the 
case here. 



Fig. 4. Residuals of Th and FPI positions with respect to a simple 
polynomial wavelength solution for HARPS orders 116 - 135 (black 
dots). The difference of the wavelength solutions between the polyno¬ 
mial model and the model with variable pixel sizes (red lines) shows the 
captured systematics. 
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Fig. 5. Effective cavity width measurement of the HARPS FPI using the 
Thorium wavelength solution. Upper panel: Measured effective cavity 
width (black dots), spline fit to red and blue HARPS CCDs (red and 
blue line, respectively). Lower panel: Residuals of the spline fit. Orders 
89 - 153 are shown here. 


This HCL wavelength solution was used as input for the FPI 
calibration to assign wavelengths to all FPI peaks. We chose 
m i = 27526 (di = 5304.24 A) which resulted in the minimum 
penetration depth variation. The model for d(k ) was a B-spline 
with 23 nodes in the blue and 13 nodes in the red HARPS regime. 
The effective cavity width measurement dk = 2 At (mi + k) and 
the according model d(k) are shown in Fig. [3 

The mean effective cavity width we measure for the HARPS 
FPI is J = 7.30020 + 0.000 15 mm. (The spike seen around 
430 nm is caused by a CCD defect in the blue detector which 
systematically influences the line position measurement in this 
region.) 

From the B-spline model and the peak number we derived 
the calibrated FPI wavelengths. To check the quality of the Tho¬ 
rium wavelength solution we plot the residuals of the calibrated 
FPI lines included in block 2 of the blue detector (echelle orders 
116 to 135) in the lower panel of Fig. [4] Due to the numerous 
more lines offered by the FPI we are now able to see systematics 
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that were hidden in the Th residuals simply due to the lack of 
lines. 

These systematics are already known from laser frequency 
comb calibration and are due to a stitching effect dWilken et al.1 
l2010t) . The 4k x 4k HARPS detector consists of two 4k x 2k 
CCDs; each is build up of 8 x 2 blocks with 512 x 1024 pixels. 

To capture the effect caused by the CCD imperfections in 
our wavelength solution, we applied our pixel size concept (see 
Sect. 3.2) and model the size of every 512th pixel as indicated 
in Eq. (|5j. From the HARPS flat fields we see that the first pixel 
block is 3 pixels smaller than the rest and ends at pixel 509 while 
the rest continues in 512 pixel intervals. The edge of the detector 
is a small block of 3 pixels which is not fitted. Furthermore, we 
assume that the last pixels in the same column have the same size 
which allows us to limit the number of free pixel parameters to 
14 per CCD. We found it necessary to increase the polynomial 
orders (8th degree in dispersion direction and 7th degree in cross 
dispersion). This is possible due to the numerous FPI lines. 

We overplot the model difference between the simple poly¬ 
nomial and the full model including inhomogeneous column 
sizes as solid red lines in the lower panel of Fig. [4] Our new 
wavelength solution handles most of the systematics a simple 
polynomial model cannot account for. Hence using the FPI along 
with a wavelength solution including unequal pixel or column 
sizes is a powerful method to characterize the detector and cali¬ 
brate the spectrograph at the same time. 

TableQ]shows the deviation of column sizes from unity, c/ plx - 
1. The smallest column we found differs by 5.8% is size which 
translates to about 46 m/s for HARPS. 

4.2. Comparison to DRS 

To check the accuracy of our calibration using the HARPS FPI 
we compare our wavelength solution to the standard calibration 
delivered by HARPS DRS (Data Reduction Software). The dif¬ 
ference between our wavelength solution and the cubic polyno¬ 
mial reverse regression of DRS is shown in the upper panel of 
Fig. [6] We observe distortions with amplitudes of up to 50 m/s 
that are repeat ing in every order. T his is very simila r to th e re¬ 
sults found in IWilken et al.1 (1201 Oh and IMolaro et al,l (1201 3h us¬ 
ing FFC data. These authors argue that the pattern originates 
from DRS not taking into account different sizes of detector pix¬ 
els. The difference between the FFC solution of IMolaro et al.1 
( 2013 1) and the DRS solution yields an rms scatter of 25 m/s. We 
confirm the distortions in the DRS solution with the use of the 
HARPS FPI and find an rms scatter between our wavelength so¬ 
lution and the DRS solution of 22 m/s. 

4.3. Comparison to LFC 

It is also interesting to compare our wavelength solution di¬ 
rectly to the FFC solution. Unfortunately, there is no night in 
the HARPS archive with FPI and FFC spectra. Hence we com¬ 
pute the difference between the FFC wavelength solution of 
IMolaro et al.1 d2013l) for November 24, 2010 and our FPI wave¬ 
length solution for March 16, 2011 and correct for 6 m/s which 
is the mean drift between both Thorium solutions. The result is 
plotted in the lower panel of Fig. [6] 

The deviations between our FPI wavelength solution and the 
FFC solution of IMolaro et al.1 d2013l) are smaller, only 8 m/s. 
Generally, no large amplitude distortions of 50 m/s are visible 
as it is the case with the ThAr solution. 



wavelength A [nm] 

Fig. 6. Upper panel: Difference between our FPI wavelength solution 
and the DRS wavelength solution (solid lines) with an rms of 22 m/s. 
Difference between FPI wavelengths obtained from the model d{k) and 
the wavelengths assigned to the FPI from DRS wavelength solution 
(points). Lower panel: Same as above but using the LFC wavelength 
solution of IMolaro et al.1 d2013l) . Orders 116 - 128 are shown here. 


We comment now on some remaining features. The peak 
around 477 nm is due to the low flux level of the LFC at this 
wavelength. 

In six pixel wide regions at the block borders (marked in cyan 
and magenta in the lower panel of Fig. [6]) the deviations can jump 
by 50 m/s. This is due to an indexing problem. We realized that 
the discontinuities (with 5jJ pixe l perio d) start at pixel 515 in 
the LFC solution of IMolaro et al.l d2013l) . We chose 509 as start 
pixel to match the columns having anomalous sensitivity in the 
flat echelle image likely due to their pixel size. 

Apart from this deviations we observe a residual trend across 
the orders which is not well understood, but we note that four 
months passed between both solutions. Around this trend differ¬ 
ences with amplitudes of about 10 m/s are seen. 

With the current data quality of the HARPS FPI spectra we 
cannot completely resolve all the fine structure that is seen with 
the LFC. IMolaro et al. ( 2013 ) fit a cubic polynomial for each 
block in each order individually, i.e. their model has more pa¬ 
rameters and is capable of handling more small scale structure. 
We used a less flexible model due to the problems with the 
HARPS FPI spectra discussed in detail in Sect. 4.4. If the data 
quality improves, we might also be able to use more detailed 
models and to characterize the detector in more detail. 

4.4. Systematic high frequency noise in the HARPS etalon 

When we calibrated the HARPS FPI the rms scatter around the 
B-spline model of cl(k) was 25 m/s, see Fig. 0 This is about 
one order of magnitude more than what we expect from photon 
noise for signal to noise ratios of 150 to 300. We found that the 
residuals are not white but exhibit high frequency variations. The 
position residuals of lines present in adjacent and overlapping 
orders are correlated indicating a wavelength dependency. Fre¬ 
quency analysis revealed two significant periods. For order 131 
we found Pj = 1.5056 ± 0.0017 A and P 2 = 0.6612 + 0.0005 A. 
Furthermore, periodic peak height variations of up to 10% are 
visible in the ID extracted FPI spectra. Frequency analysis on 
the derived FPI amplitudes yields two significant periods match¬ 
ing those found in the position residuals. We presume that an 
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Table 1. Deviation of column sizes, a pix , from unity in percent for F1ARPS. Estimated errors are typically 0.3% of a pixel (~ 2.5 m/s). 


Column number i 

509 

1021 

1533 

2045 

2557 

3069 

3581 

Block row 1 (o = 135 - 161) 

-4.3 

1.1 

-1.1 

-3.5 

0.5 

-4.5 

1.3 

Block row 2 (o = 116- 135) 

-2.2 

-0.8 

-1.4 

-1.0 

-1.5 

-5.8 

-0.5 

Block row 3 (o — 100 - 114) 

-5.3 

2.4 

-3.1 

-3.4 

-0.5 

-4.6 

1.6 

Block row 4 (o — 89 - 99) 

-2.8 

-1.0 

-5.0 

-0.4 

0.2 

-5.8 

-0.2 


additional low finess interference signals, caused by plane par¬ 
allel optical elements in the path of the FPI like filters or by the 
thickness of the FPI mirrors itself, superimposes on the HARPS 
FPI (tWildi et al.l (1201 Ot) . Pepe, private communication). The pe¬ 
riods would correspond to FPIs with cavity widths in the order 
of millimeters. Further we find that the periods decrease towards 
blue and increase towards red orders, indicating that the effective 
cavity width of the optical element causing the additional FPI 
effect changes with wavelength; this can be caused by a wave¬ 
length dependent refractive index n. We presume that the slopes 
due to the flux variations shift the photo-center of the FPI lines 
resulting in systematics in the line position measurement. 

Without exact knowledge about the cause of this effect, our 
attempts to fit and correct the flux variations present in the spec¬ 
tra gave unsatisfying results. Therefore the FPI wavelengths are 
not known precise enough to resolve the fine structure that is 
seen in the LFC. Our method involving the FPI in the wave¬ 
length calibration for HARPS is thus limited by this systematic 
high frequency noise. Our laboratory experiment in Sect. 2.2 is 
limited by photon noise (8 m/s). With this data quality we would 
be able to resolve more details in the wavelength solution sys¬ 
tematics of HARPS. For future projects it is advisable to avoid 
optical elements in the light path of the FPI that might cause 
additional interference patterns and to characterize the FPI in 
advance with other instruments. 


5. Influence of distortions in the wavelength 
solution on precise RV measurements 

To estimate the influence of distortions in the wavelength solu¬ 
tion on precise RV measurements with high resolution echelle 
spectrographs we simulate observations taken over the period 
of one year. The simulated spectrograph has 70 spectral orders 
each covering 57 A, a delta peak instrumental profile and a de¬ 
tector with 4096 homogeneous pixels. The observed star is per¬ 
fectly quiet and has a zero RV signal. Its spectrum consists of 
randomly distributed delta peaks (100 lines per order). Further¬ 
more, the star lies in the ecliptic resulting in a 30 000 m/s RV 
signal from the barycentric motion of Earth. This oversimplified 
approach allows us to single out the effect of the wavelength so¬ 
lution only on our observations without additional noise from 
any other sources. 

In total we simulate 29 observations over one year. We ini¬ 
tialize the stellar spectrum by computing random wavelengths 
such that 100 delta peaks fall on each of the 70 orders of the 
spectrograph. At time stamp 0 the barycentric motion of Earth 
towards the star is zero; hence this spectrum represents the rest 
frame. The spectra for all observation time stamps are computed 
applying only the Doppler shift corresponding to the barycen¬ 
tric motion of Earth around the Sun (~ 30 000 m/s). We ignore 
the effect of Earths rotation because it is small compared to the 
motion of Earth around the Sun (~ 500 m/s). Spectra are taken 
every 12.175 days. We use the inverted wavelength solution of 



Fig. 7. Simulation of the RV signal due to a disturbed wavelength so¬ 
lution. RV signal for single orders (gray dots), mean RV of all orders 
(black squares), sinusoidal fit to the mean RV (red line) 


the spectrograph to place the stellar lines on the detector and to 
obtain our simulated observations. 

If the true wavelength solution of the simulated spectrograph 
is used to convert the pixel value back to wavelength and the 
barycentric motion of Earth is corrected, we obtain RVs of 0 m/s 
for all stellar lines within machine precision. 

Next we disturb the wavelength solution of the spectrograph. 
We want to test the effect of small systematics in the wavelength 
solution on precise RV measurements. Hence, the difference be¬ 
tween our Thorium wavelength solution and our calibrated FPI 
wavelength solution in HARPS order 118 (Sect. 4.1, red line in 
Fig ® is added to the true wavelength solution of the simulated 
spectrograph. 

Now we use the disturbed wavelength solution, convert the 
pixel values into wavelengths, correct for the barycentric motion 
of Earth and compute the RVs for all lines. We average the RV 
of single lines to compute the RV order by order for each obser¬ 
vation (gray dots in Fig.®. The RV per time stamp is the mean 
of the individual order RVs (black squares in Fig.®. 

We observe the following effects in Fig.® The RVs derived 
for individual orders (gray dots) show a scatter of about 1 m/s. 
This effect is a result of the random distribution of lines on the 
CCD and the systematics in the wavelength solution. 

Depending on the location of a stellar line on the CCD it 
will give a different systematically wrong RV because the sys¬ 
tematic error in the wavelength solution is different for different 
pixels. We can derive the expected RV signal S v resulting from 
the wavelength solution systematics djx) by inserting the wave¬ 
length solution distortion in the Doppler equation: 

[/1(,Y/) + 6Axd] - [d(,Yo) + diCro)] _ v + 6 V 
Mx 0 ) + 6 A (x 0 ) c 

where xq and x; correspond to the line positions on the CCD at 
time stamps to = 0 (rest frame) and r, respectively. The barycen¬ 
tric motion of Earth shifts lines on scales of a few 10 pixels. On 
this scale the change in d,\ is small and we can approximate the 
difference S^x i) - S,i(xo) with the derivative of the wavelength 
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solution distortion multiplied by the pixel shift dd,i(x)/d.\ ■ Ax. 
The pixel shift Ax depends on the barycentric velocity of Earth 
RVe and the pixel scale in terms of RV k. For HARPS one pixel 
covers approximately 800 m/s/pix. Hence a barycentric velocity 
of RVe = 30000 m/s shifts the lines by 37.5 pixels on the de¬ 
tector. Finally the small term c)j(xo) in the denominator can be 
neglected and velocity signal of a single line <5, due to the distor¬ 
tion in the wavelength solution djx) can be written as 


S v 


c d6,\(x) 1 Rv ^ 
A(x 0 ) dx k 


( 8 ) 


Averaging over all the different 6 V in one order results in a sys¬ 
tematic RV that depends on the random distribution of stellar 
lines on the CCD. As the distortion is the same for all orders, 
the simulation results presented in Fig. [7] are equivalent to re¬ 
peating the test 70 times for one order with different random sets 
of lines. Hence we observe scatter between the RVs of different 
orders only caused by the line distribution on the detector. This 
scatter is mainly caused by lines experiencing different slopes in 
&5,\(x)/Ax rather than by single lines moving across the jumps in 
the wavelength solution. 

The slope dd,i(x)/dx is, however negative almost everywhere 
on the CCD. Hence many lines experience a systematic RV shifts 
in the same direction. Therefore, the mean RV derived from the 
individual orders shows a low amplitude RV signal with a period 
of one year. The amplitude of this signal is about 0.5 m/s. Ad¬ 
vancing to cm/s precision we cannot afford systematic errors in 
the order of 0.5 m/s to 1 m/s. Hence, accurate wavelength cali¬ 
bration is critical for reaching the aims with the next generation 
of instruments. 


- Motivated by the grating equation we parametrize our model 
as x(A, o ) = poly(o/l, o) which is a better representation of 
the problem than poly(/l, o). 

- We account for inhomogeneous detector pixel sizes a; and 
transform the observed pixel positions X to focal plane coor¬ 
dinates x(X, a ;) before performing direct regression. 

Using this concept for the HARPS FPI we are able to char¬ 
acterize the detector by taking into account pixel stitching and 
compute the wavelength solution simultaneously. 

The standard DRS HCF calibration deviates from the FFC 
solution with an rms scatter of 25 m/s. Direct comparison of our 
FPI wavelength solution with the FFC yields only small differ¬ 
ences of rms = 8 m/s. Hence we show that the calibrated FPI is 
a suitable calibration source for high resolution echelle spectro¬ 
graphs. 

Finally, we performed simulations to derive the fundamental 
limits an imperfect wavelength solution imposes on RV mea¬ 
surements and find that even small distortions in the wavelength 
solution lead to systematics signals and additional noise, both 
in the order of 1 m/s. We conclude that wavelength calibrators, 
such as FFCs or FPIs, are necessary that provide a much denser 
line grid than HCFs do. We suggest that a calibrated FPI can 
be considered as economical alternative to FFCs to achieve very 
high RV precision for current and future instruments. 
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6. Summary and conclusions 

HCFs have been calibration sources for echelle spectrographs 
since decades but advancing to cm/s RV precision introduces the 
need for more suitable standards. Aside from the FFC, FPIs also 
offer a dense comb pattern of lines that can be used for calibra¬ 
tion. 

Because the cavity width is not known accurately enough and 
possibly drifting over time, the peak wavelengths are uncertain 
and FPIs cannot be used as a standalone calibrator. Therefore 
we developed a method to anchor the dense etalon lines with the 
absolute HCFs to combine their precision and accuracy. Cen¬ 
tral to our method is that we model the cavity width (includ¬ 
ing its global wavelength dependency) in the calibration pro¬ 
cedure. We exploit the assumption that the wavelength depen¬ 
dency of the cavity width is smooth which enhances the quality 
of the wavelength solution on medium to small scales, while the 
global accuracy is still determined by the HCF. The assumption 
of smooth cavity width was verified with FTS measurements on 
CARMENES FPI. 

Along with our model for the FPI, we described our model 
for the wavelength solution of the spectrograph. Our four main 
points for the calibration are: 

- We suggest to perform direct regression with x(A) rather than 
reverse regression A(x) because the independent quantity is 
the wavelength A and not the pixel position x of the spectral 
line on the detector. 

- We use a 2D polynomial as model in the direct regression. 
This couples orders, reduces the number of parameters and 
enhances the robustness of the fit. 
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